ULTRASPHERICAL TYPE GENERATING FUNCTIONS FOR 
ORTHOGONAL POLYNOMIALS 



Abstract. We characterize, under some technical assumptions and up to a 
conjecture, probability distributions of finite all order moments with ultras- 
pherical type generating functions for orthogonal polynomials. Our method is 
based on differential equations and the obtained measures are particular Beta 
distributions. We actually recover the free Meixner family of probability dis- 
tributions so that our method gives a new approach to the characterization of 
free Meixner distributions. 

1. Motivation: Meixner families 

There is a one to one correspondance between probability distributions on the 
real line and polynomials of a one variable satisfying a three-terms recurrence rela- 
tion subject to some positivity conditions ([9]). That is why in most of the cases, if 
not all, one tries to characterize probability distributions using generating functions 
for orthogonal polynomials. Among the famous generating functions are the ones 
of exponential type, that is if fj, is a probability distribution with a finite exponential 
moment in a neighborhood of zero 



c fi(dx) < 



then 



(1) i>(z,x):=YP n { 



n>0 



e xH{z) 
E ( e Xi?(*))> 



where H is analytic around z — such that H(0) — 0,H'(Q) = 1, X is a random 
variable in some probability space (fi, & , P) with law \i = PoA -1 and (P n )n>o is the 
set of orthogonal polynomials with respect to \x. Up to translations and dilations, 
there are six probability distributions which form the so-called Meixner family 
referring to its first appearance with J. Meixner ([2]). It consists of Gaussian, 
Poisson, Gamma, negative binomial, Meixner and binomial distributions. This 
family appeared many times under differents guises ([E], [13], PQ, [15], [TT]). 
Another well known example was first suggested and studied in [2] and is given by 
a Cauchy-Stieltjes type kernel. Namely, if fi is a probability distribution of finite 
all order moments, then 

(2) *(Z, X) := J2 P n(x> n = , , sri , T 

^ Hz)[f{z)-x] 
where u and z i— » zf(z) are analytic functions around zero such that 

lim U ^ = lim zf(z) = 1. 



Key words and phrases. Generating functions, ultraspherical type, orthogonal polynomials, 
Jacobi-Szego parameters. 
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This family, known as the free Meixner family due to its intimate relation to free 
probability theory, covers six compactly-supported probability measures too. We 
refer the reader to [3], [5], [5], [12j for more characterizations and more interpreta- 
tions. The natural g-deformation that interpolates the forementioned families for 
arbitrary \q\ < 1 was defined and studied in [3] and is up to affine transformations 
the so-called Al-Salam and Chihara family of orthogonal polynomials ([I]). Their 
generating functions is given by an infinite product and is somehow similar to the 
^-exponential function. Another characterization of the last family was recently 
given in [TJ. 

After this sketchy overview, we suggest another type of generating functions which 
may be viewed as a generalization of the free Meixner family. It is inspired from 
the case of Gegenbauer or ultraspherical polynomials for which ([9]) 

(3) <£ 2 »®?-Cfa)z» = 1 M < 1,A > 0, 

' (1 — 2zx + z z ) A 



n>Q 



where (A)„ = (A + n — 1) • • • (A + 1)A and for complex z such that the RHS makes 
sense and the series in the LHS converges. We adopted here the monic normalization 
for (C£)n an d henceforth all the polynomials are monic so that they satisfy the 
normalized recurrence relation 

(4) xP n (x) = P n +i(x) +a n P n {x) +uj n P n - 1 (x), n > 0,P_i := 0,w = 1- 

The sequences (a n ) n >0y {w n )n>o are known as the Jacobi-Szego parameters and 
u> n > for all n unless \i is has a finite support ([§])■ Moreover, we shall always use 
these notations for the different families of orthogonal polynomials we shall cross 
through this paper. 

It is then natural to adress the problem of characterizing probability measures of 
finite all order moments, say fj,\, such that 

(5) : =g o ^^"= u xm X)-*r A> °' 

valid for x £ supp(/x,\) and z belongs to a complex open region S near z = cut 
from z = along the negative real axis where u\, fx are analytic with 

(*) lim zh(z) = 1, lim = 1, 9f(/(*)) /0,zeS. 

zes 

By the last assumption, (f(z) — a;) A is well defined for all x € supp(/z), z £ S and 
A > (the principal determination of the Logarithm is adopted). Moreover, the 
above limiting conditions imply that t^\(z, x) tends to 1 as z tends to in S for all 
x £ supp(/ix). We shall say that ip\ is a generating function for orthogonal polyno- 
mials of ultraspherical-type referring to ultraspherical polynomials. Without loss 
of generality, we may assume that fix is standard, that is, has a zero mean and a 
unit variance. Equivalently, if (a A )„>o, (w^) n >o denote the Jacobi-Szego parame- 
ters of fi\, then one has ckq = 0, to* = 1. Our strategy is based on the following 
general claim that was stated without proof in ,6j and proved below for the reader's 
convenience: 

Claim: to a given generating function for orthogonal polynomials (z, x) i— > ip(z, x) 
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associated with a (standard) probability measure fi satisfying some integrability con- 
ditions (to be precise later), the measures {V z } defined by 



V z (dx) := ip(z, x)/j,(dx) 

are probability measures such that the mean and the variance of P 2 are polynomials 
in z of degree 1, 2 respectively. {P 2 } is then referred to as the ^-family of fi with 
an at most quadratic variance, referring to both the exponential and the Cauchy- 
Stieltjes families ([8], [15]). When ip is handable enough so that one can perform 
computations of the first and of the second moments of F z independently from the 
infinite series, one recovers two equations that may be used to solve the problem of 
characterization of probability measures whose generating function for orthogonal 
polynomials is given by ip (or of ?/>-type). In the case of the Meixner and the free 
Meixner families, this was noticed in [5J. In the case in hands, if the assumptions 
in (*) are valid for z € S together with the assumption (**) (see below) , we obtain 

Proposition 1.1. 

(1) The function f\ satisfies for z£S 

(6) Q 2 (z)f' x (z) = f 2 x {z) - Qi(z)f x (z) + R^z) 

where Q2, Ri are polynomials of degree 2 while Q\ is a polynomial of degree 
1. Moreover the coefficients of these polynomials depend only on A, a^,^- 

(2) The function u\ is related to f\ by 

«A(*) _ x 1 - f'xW 



X- 



y(z) E x {z) 



u\(z) fx(z)-\z 
Once we did, we show that if 

(7) gx(z) := h(z) - 

where E\ is assumed to be a polynomial, then deg(E\) < 2 and this follows from 
the fact that Q2,Qi, R\ are polynomials (terminating series). Next, we investigate 
under the last assumption the case of symmetric measures. We show that there 
exist two families of probability measures corresponding to (C*)„ for A > and 
(C„ *)« for A > 1/2, A 7^ 1. We warn the reader to the fact that, though these 
two families differ from each other by a parameter's translation, their generating 
functions given by ([5]) are totally different since depends on A and is fixed 
for both families. Under the same assumption, there is only one family of non 
symmetric probability measures corresponding to shifted monic Jacobi polynomials 
p A-i/2,A-3/2 jP A-3/2,A-i/2 for x > y 2> X^ 1. The discard of the value A = 1 is 
needed for the computations since we need to remove factors like 1 — A, 1 — A 2 . 
Thus, one deals with this case separately and recovers the free Meixner family for 
which deg(i5>,) < 1 too. 

Problems: we do not know if there exists a solution f\ for which E\ is an entire 
infinite series. Note that such a solution does not exist when A = 1. However, 
we already know that the free Meixner family covers six families of probability 
distributions ([4]) while there are three families for A 7^ 1 when E\ is a polynomial. 
Is there any intuitive explanation to this difference between both cases or to the 
degeneracy of the case A = 1? 
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n>0 n>0 

ix)z n = aiu>i = a%z. 



2. Validity and Proof of the claim 

Write tfi as 

ip(z,x) = a n P n (x)z n 

n>0 

for some fixed sequence (a n ) n ,x G supp(/i) and z in a suitable complex domain D 
near z = so that the infinite series converge. The integrability conditions we need 
for the claim to be valid are the finitcness of all order moments of /z and 

f a n {x l P n {x))z n ^{dx) = J2 a n f x* P n {x)fi{dx)z n , i G {0, 1, 2}, (**) 

J n>0 n>0 

for z G D. In fact, for i = 0, the orthogonality of P n shows that P z is a probability 
measure for all z G D (remember that Pq = 1) and together with ao — Q, oji = 1 
imply 

a n / P n +i(x)n(dx)z n = 0, n > 0, a n a„ / P n (x)fi(dx)z n = a Q a = 0, 

ra>0 

^a„w„ / P n ^i(x)ix(dx 
Thus, one gets for i = 1 after using (l) 

(8) /" ^ a„(xP n (x))zXcfe) = a lZ = / xP,(dx). 

n>0 

For i= 2, one uses twice (j4} to get 

(9) / ^ a„(x 2 P„(x))z"^(dx) = a 2 w 2 2 2 + a x a x z + 1 = / £ 2 P 2 (efe) 

n>0 ^ 

and the claim is proved. 

Remark 2.1. In the case in hands, if [i\ is compactly supported, then the Jacobi- 
Szego parameters are bounded thereby one can exchange the infinite sum and integral 
signs. Indeed, by Cauchy-Schwarz inequality 

E^T / WPn(x)\»x{dx)\z\ n < (f |x| 2 v(dx)) 1/2 ^^i|p„nizr, 

ri>0 ' J \J J „> 

for i G {0,1,2}. Moreover, ||P r i|| 2 = ■ ■ •k'n-i < c " for some c > so that 
Fubini's Theorem applies for \z\ < 1/s/c. As the reader can see, the exchange of 
the order of integration depends on the sequence (a n ) n and the growth conditions 
satisfied by /i. As a matter of fact, if («„)„ is fixed, they solely depend on [i (or in 
\\Pn\\). 

3. Proof of Proposition II. II 

3.1. First and second moments. On the one hand, the integration of both sides 
of ([5|) with respect to fi\ gives 
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On the other hand, one gets from (JSJ), © and a n = (A) n /n! 



m 2 (z) := / x tjjx(z, x)fi(dx) — lo 2 z + Aaqz + 1. 

Then, using the elematary operation x = (x — f(z)) + f(z), it follows that 

m " (z) = f{z) ~ UA4(2) ^ i {f{z)-xY-^ x[dx) - 

Differentiating with respect to z G S under the integral sign0 defining Ux,i, one 
gets (1 — X)f (z)u\(z) = (uxi)'(z). Thus the RHS of (z) transforms to: 

(10) <M = A . 1 -/aW 



u\(z) fx(z)-Xz 
which can be written as 

(11) (u x (z)[f x (z) - Xz)})' = (1 - X)ux(z)fx(z). 
For the second moment, use x 2 — x(x — f(z)) + xf(z) to get 

(12) ml(z) = Xzfx(z) - ^ I {h(z) X _ x) x^x(dz). 
Using 



—Hxidx^j = (l-X)f'x(z) jj-jJj—-^Lix{dx) = X(l-X)zux{z)f' x {z) 



(fT2| is rewritten as 

(13) {[Xzfx(z) - m x 2 (z)]ux(z))' = A(l - X)zux(z)f x (z). 

3.2. A non linear differential equation. By the virtue of (fTT|) . (fT3|) implies that 
([AzA(z) - ro£(*)]u A (*))' = Az( WA (z)[/ A (z) - A*)])' 

which gives 

[Xzfx(z) - m x 2 {z)]u' x (z) + [Xfx(z) + Az/ A (z) - (m$)'(z)]ux(z) 
= Xz[f x (z) - Xz]u' x (z) + Xz[f x (z) - X)]u x (z), 

therefore 

[X 2 z 2 - m$(z)]u' x (z) = [K)'(z) - Xfx(z) - X 2 z]ux{z). 
If Xz — rn^fz) 7^ 0, one gets after the comparison of the last equality to (TT0|) 

(m$y(z)-Xfx(z)-X 2 z } l-f x (z) 
X 2 z 2 — m\(z) f\(z) — Xz 

which shows after elemantary computations that fx satisfies the following non linear 
first order differential equation: 

(14) Q2(z)f' x (z) = fl(z) - Qi(z)fx(z) + Ri{z) 



This is justified by the analyticity of fx in S and general properties of generalized Cauchy- 
Stieltjes transforms, see |17j and references therein. 
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where 



Q 2 (z) = x 



A - h±lu* 

2 2 

A. 



z 2 - \ot\z - 1, 



Q x (z) = (\+l)cj$z + al, 
Rl {z) = ^±^ 2 V-1- 
Setting gx(z) := / A (z) - [Qi(z)/2], |14| transforms to 

(15) Q2(z)g' x (z) = g 2 x (z) + Q 2 (z) 
where 

Q 2 (z) = Ri(z) - \[Qi{z)f - ^±l^Q 2 (z) 

. \ . A 2 — 1 A 2 A 2 — 1 % \ (A + 1)lj 2 a (a A ) 2 
= [(A + l)c A - 2A]— ^a,*z 2 + — a^z + i ^ - 1 - 

Finally, once g\ is given, one deduces f\ by adding Qi/2 then use (jTUJ) to derive 

U\. 

4. Some solutions of © 
From now on, we shall look for solutions of ^ of the form 

gx{z) : = EM, Ex (Q) = l 

for a second degree polynomial E\. In fact, since z <— > zgx(z) is analytic around 
zero, one may always assume that g\{z) has the above form for an entire function 
E\. But if E\ is a polynomial of degree > 3, then all the terms of degree > 3 will 
vanish only by equating both sides of (fl~5|) . For instance, let 

E\(z) = a z 3 + aiz 2 + a 2 z + a 3 

and write (| 1 5|) as 

(16) Q 2 {z)[zE' x (z) - E x {z)] - El(z) = z 2 Q 2 (z). 

Then by equating terms of degree 6 is this equation, one easily gets ao = so that E\ 
has degree 2. For E\ a polynomial of degree 4, start with equating terms of degree 
8 and so on. However, this way of thinking fails or rather become cumbersome 
when E\ is an entire function and the existence of such a solution is open. 

4.1. A new approach to the Free Meixner family. Recall that the free Meixncr 
family corresponds to A = 1 and that it covers six compactly-supported probability 
distributions given by their Jacobi-Szego parameters ([!]) 

ct\ = a, a el,n> 1, = (1 + &),&> -l,n > 2, 

where we used the fact that \x\ has a mean zero (a\ — 0) and a unit variance 
(io\ = 1). Moreover, one has ([5]) 

, , , l + az + (l + b)z 2 . . (a/2)z + l a 1 

z z 2 z 

But, Q 2 reduces to a constant for A = 1 so that (fT5|) transforms to 

[(1 - ^)z 2 a\z - l]g' x (z) = g 2 (z) + (c A 1) - (a A ) 2 /4. 



It is then an easy exercice to check that g\ satisfies (fT5|) which reads in this case 

(17) - [bz 2 + az + l]g[{z) = g\{z) + b - a 2 /A. 

We can even prove that g\ as written above is the unique solution of the last 
differential equation subject to the condition zg\(z) — > 1 when z — ► 0. In fact, 
writing g\(z) — h\(z) + l/z is some punctured neighborhood of zero where hi is 
analytic around zero, simple manipulations show that hi satisfies 

-[bz 2 + az + l}h[(z) = h\{z) -^ + - z (hi(z) - ~) . 

Taking the limit as z — > 0, one has from the singularity at z — in the RHS that 
h\{Q) = a/2. Thus, one gets 



[bz 2 + az + 1] ^ nc n z n 1 = . 



CnZ 

n>l n>l 



^ c n z n + a 



n>l 



for some sequence (c„) n >i, which makes sense for z = therefore c\ = 0. Removing 
z from both sides of the obtained equation then setting z = will give c 2 = 0, 
removing z 2 and taking z — gives C3 = and so on. As a result, h\(z) = a/2 
and our method gives a new (geometrical) approach to the characterization of free 
Meixner distributions. 

Remark 4.1. When A 7^ 1, auxiliary terms show up and h\ satisfies 

Q 2 (z)h' x (z) = h\{z) - + I (hx{z) - X^j + ^(2 - W *) 

which shows that h\(0) = Aa^/2 w/izZe h'^(0) ^ 0, /i"(0) 7^ in general. 

5. Symmetric measures: ultraspherical polynomials 

In the sequel, we shall focus on the case = for all n. This is equivalent to 
the fact that fi\ is symmetric, that is the image of fi\ by the map x 1— > —a; is still 
/ia- In this case, one gets by taking a\ = 

Q a (z) = ^[2A-(A + lK]z 2 -l, 

Q 2 (z) = [(A + lK-2A]^^ + ^M-l. 

Writing E\{z) — oqz 2 + aiz + a 2 and equating both sides in (|15[) . one gets: 

a 2 = 1, 
ai = 0, 

_3ao-^[2A-(A + lK] = - 1, 

A A 2 — 1 

-a 2 + a -[2A- (A + 1V 2 A ] = [(A + 1)^ 2 — 2A] — - — u)%. 

The third equation gives 

(1-A 2 )(2-^) 

6 — T' 

Hence, it remains to check when the above ao satisfies the fourth equation. Since 
the case A = 1 is known, we assume A ^ 1 so that one removes the term (1 — A 2 ) 
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in the above equalities. Substituting ao in the fourth equation, one sees that uj^ 
satisfies 

-(A + 1)(A + 2)(c^) 2 + (4A 2 + 6A - + (1 - 4A 2 ) = 0. 

What is quite interesting and even surprising, that though this polynomial looks 
complicated, its descriminant is equal 9 so that there are two solutions given by 

2A+1 , 2A-1 



A+2 ^ A+l 

where for the second value, we consider A > 1/2 in order to avoid finitely-supported 
probability measures and signed measures. As a result, 

1 - A 2 1 - A 2 1 - A 



a — 7771 — ; TTTi a — 



Thus 



and from fl" 



Finally 



2(A + 2)' 2(A+1) 



f f \ 1 + A J- 1 f t \ X JL 1 



u' x {z) _ A u' x (z) _ A z 2 + l-(A/2)z 2 
u x {z) ~ z' u x {z) ~~ z(l-(A/2)z 2 ) ' 



z x 

u x {z) = z\ A > 0, A ^ 1, u x {z) = - _ pr^ .A > 1/2, A ^ 1, 



for z E S. Note that 5 is easily described: in fact /a is not real outside the real 
line and the circle \z\ < 2/(1 + A) or \z\ < 2/A respectively. Moreover the fi\ is 
compactly-supported as we shall see below, so that (**) is satisfied in a ball centered 
at the origin (see remark |2~T]) . 

5.1. Ultraspherical polynomials: symmetric Beta distributions. The value 
i corresponds to the ultraspherical polynomials. However, in order to fit into 

our setting, one has to consider the monic Gegenbauer polynomials, say C A , which 
are orthogonal with respect to the standard Beta distribution 

c A (l-2; 2 /[2(l + A)]) A - 1/2 cfa, x € [±\/2(l + A)] 
for some normalizing constant c\. They are given by 



= (\/2(TTA)rC„ 



V2(l + A) / 
Now, it is easy to see from ([3]) that 



n 

n>0 



1 



(l-zir+(l + A)z 2 /2) A 



l + (l + A)z 2 /2 



u x (z)(f x (z) - x)*- 



For w^, i>\ is written as: 

ip\{z,x) 



1 - (A/2)z 2 



1 - (A/2)z 2 



z x (\z/2 + l/z-x) x (Az 2 /2 + l-zx) A 

and we claim that P A = C n for all n and all A > 1/2, A 7^ 1. In fact, 
(A) n/5 A-l, 



n>0 



n>0 



A + n-l(A-l)» - A-i 
A-l n! °" lj 



i V — —C X ~ 1 (x)z n+X - 1 



n>0 



n! 



1 



(A - l)z A ~ 2 

1 - (A/2)z 2 
(l-zx + Az 2 /2) ; 



l-zx + Az 2 /2 



A-l 



as the reader may easily check. 



6. NON-SYMMETRIC PROBABILITY MEASURES: JACOBI POLYNOMIALS 

Henceforth, we suppose that a x 7^ 0, A 7^ 1 and we will show that there is only 
one family of probability measures subject to 



9\{z) = 

Then, we get the following equations 



a 2; + a\z + a 2 



«2 

ai 



-3a --[2A-(A + l)^]-a 2 = 



(A + 1 K) 



A\2 



- 1 - 



-aoftiA — 2a fti 



A 2 -l 



a A w A , 



A A 2 — 1 

-a 2 + a -[2A-(A + lV 2 A ] = [(A + 1)^ - 2A] — - — uj^. 

From the second, third and fourth equations, it follows that 



a 



1-A 2 
6 



1-A 2 
4A 



Actually, this gives a constraint on A,a A ,w A; 
(18) <>» 2 



2 A 



A + ^ ) ^2- 



Substituting ao by (1— A )w2/(4A) and removing (1 — A ), the fifth equation becomes 

~ )2 + T [2A " (A + = [2A " (A + ^IT 



In the non degenerate case ^ 0, 



4A 3 



" 2 - 2A3 + 3A 2 -1' 
But —1 is a double root of the polynomial in the denominator so that 

2A 3 



2 (A + l) 2 (A-l/2)' 
which is positive for A > 1/2. Finally, one deduces from (fl8|) that 



(^) 2 = 2 

a = = 
It follows that 

h(z) = 



(2A + 3)A 2 



- 2 



(A + l)2(A-l/2) 
(l-A)A 2 



.(A + l) 2 (A-l/2) 

(1-A 2 )A 2 

2(A + 1) 2 (A - 1/2) ~ (A + 1)(2A - 1) ' 

QqZ 2 + (l\Z + 0,2 (1 + A)o>2 z + a l 



>0, 



z 

1 - A 
2A 
_>? 
2A- 1 



1 



z 2 ± 



1 + A 
2 
1 



V2A- 1 



A + l 



z + 1 



and 



u' x {z) _ A 



1 



(A -I) 2 



2A- 1 

The descriminant of the polynomial 



Aa-A) z2± 

2A - 1 
1 



V2A- 1 



=z + 1 



is easily seen to be: 



A(l~A) r2± 
2A- 1 V2A - 1 

1 4A(1-A) 



z + 1 



2A - 1 



2A- 1 



2A - 1 > 0. 



It follows that, when a* > 0, the roots are given by 



Writing 



1 - 



(A - I) 2 9 



2A- 1 
one gets 

u' x (z) _ A-l 
As a result 



V2A- 1 
A : 



(A-l) 2 
2A- 1 



22 



V2A- 1 
1-A ' 



z + 



V2X-1 



1-A 



V2A- 1 



1-A 



V2A- 1 



A-l 



z + 



V2A - 1 
A 



A 



z z + V2A- 1/A 



V2A- 1 



A z + v/2A - 1/A 
and the generating function is written as 

V2A- 1' 



(19) t/> x (z,x) 



A 

V2A - 1 



z + 



A 



1 — z \ x 



A 2 



V2A -I J 2A-1 
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In the case a* < 0, similar computations yield 



u\(z) = - 



V2A- 1 



A z - V2A- 1/A 



and 



A 



V2A- 1 



V2A- 1 
A 



1 — z [ x 



1 



A 2 



V2A- 1/ 2A-1 



6.1. Orthogonality measures: special Jacobi polynomials. Wc will show 
that P£ is a shifted monic Jacobi polynomial with parameters depending on A. To 
proceed, recall that ([10]) the monic Jacobi polynomials are orthogonal with 
respect to the Beta distribution with density function given by 

C a j(l-x) a {l + x) fi l [ -. ltl] (x), a,(3 > -1, 

for some normalizing constant c a ,p and that the non monic Jacobi polynomials 
P"'P are related to p" ,/3 as 

n y ' 2 n n\ n K ' (a + f3 + l) n 2 n n\ n y ' 



We will show that 

>A 



2A 



when a i > and 



Pn(x) 



V2A- 1 



2A 



A-1/2.A-3/2 



/ V2A- lx- 1 



2A 



V2A- 1 



P 



A-3/2,A-l/ 



A-i/2 /V2A- lar + l 



2A 



when cti < 0. Before proceeding, note that both cases are related using P"'@(x) 



pA-3/2,A-l/2 ^ V2A- lx + 1 ^ = ( -_ 1 ^ npA -l/2,A-3/2 ^ \/2A - l(-x) - 1 ^ 

so that their generating functions are the same up to the transformation (z,x) 
(— z, — x). Moreover the orthogonality measures are given by 



H\{dx) 



H\{dx) 



ca 1 



4(1 



V2A- lx- 1 
2A 



A-l/2 



V2A - lx - 1 
2A 



A-3/2 



dx, 



y2A^Tx + l\ A_1/2 / V / 2A _ Tx + 1^ A_3/2 



2A 



1 



2A 



e?x, 



for some normalizing constants c\ , c' A and for 

1 - 2A 1 + 2A 

X £ 



x £ 



V2A- 1 V2A- 1_ 
1 + 2A 2A - 1 



V2A- l' y/2X- 1 
li 



respectively. 

Now, we proceed to the proof of our claim and we consider the case a* > 0. To 
this end, we need ([10]) 

{l + t)a+p+1 2 iy 2 2 > (i +t )2 ) (/?+!)„ » [y) 

= ST ( a + ^+ 1 hn a, 0( , ( t 

^ G9+l)„n! Pn [V> \2 

for \t\ < 1, < 1, where 2-F1 is the Gauss hypergeometric function ([9]). Sub- 
stituting (a,/3) by (A - 1/2, A - 3/2), then (a + (3 + l)/2 = A - 1/2 = /? + 1 so 
that 

/q+g+1 q+g+2 2 (y + l)^ / 2(^ + 1)^ / 2(y+l) t\- A 

where we used that iF (A,y) = (1 - y)~ x for |y| < 1 ([9]). Thus 

1 / <* + /?+! Q + ^ + 2 2(y+l)A _ 1 + t 

(l + t) Q +/3+l 2 1 ^ 2 ^ 2 ; (l + t )2 J ~ [1 + < 2_ 2 ^]A- 

Now use the Gauss duplication formula ([5]) 

V^T(2a) = 2 2a - 1 r(a)r(a + 1/2), a > 0, 



to see that 



As a result, 



(a + P + 1) 2 „ _ (2A - l) 2n _ o2n 



(/?+!)« (A-l/2)„ 



= 2^"(A) r 



1 + t 

n>0 1 aJ 



It finally remains to substitute in the last equality 

V2A - lx - 1 A 

y = 2A ' t = ^W=i z 

for small z to see that it is nothing but ([19]) and the claim follows. 
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